arXiv:1506.02025v3 [cs.CV] 4 Feb 2016 


Spatial Transformer Networks 


Max Jaderberg Karen Simonyan Andrew Zisserman Koray Kavukcuoglu 

Google DeepMind, London, UK 

{jaderberg^ simonyan^ zisserman,korayk}@google.com 


Abstract 


Convolutional Neural Networks define an exceptionally powerful class of models, 
but are still limited by the lack of ability to be spatially invariant to the input data 
in a computationally and parameter efficient manner. In this work we introduce a 
new learnable module, the Spatial Transformer, which explicitly allows the spa¬ 
tial manipulation of data within the network. This differentiable module can be 
inserted into existing convolutional architectures, giving neural networks the abil¬ 
ity to actively spatially transform feature maps, conditional on the feature map 
itself, without any extra training supervision or modification to the optimisation 
process. We show that the use of spatial transformers results in models which 
learn invariance to translation, scale, rotation and more generic warping, result¬ 
ing in state-of-the-art performance on several benchmarks, and for a number of 
classes of transformations. 


1 Introduction 

Over recent years, the landscape of computer vision has been drastically altered and pushed forward 
through the adoption of a fast, scalable, end-to-end learning framework, the Convolutional Neural 
Network (CNN) ED. Though not a recent invention, we now see a cornucopia of CNN-based 
models achieving state-of-the-art results in classification |[T3[28l[35l, localisation EDEl. semantic 
segmentation ED, and action recognition |[T2|[32l tasks, amongst others. 

A desirable property of a system which is able to reason about images is to disentangle object 
pose and part deformation from texture and shape. The introduction of local max-pooling layers in 
CNNs has helped to satisfy this property by allowing a network to be somewhat spatially invariant 
to the position of features. However, due to the typically small spatial support for max-pooling 
{e.g. 2x2 pixels) this spatial invariance is only realised over a deep hierarchy of max-pooling and 
convolutions, and the intermediate feature maps (convolutional layer activations) in a CNN are not 
actually invariant to large transformations of the input data |[6l[22l. This limitation of CNNs is due 
to having only a limited, pre-defined pooling mechanism for dealing with variations in the spatial 
arrangement of data. 

In this work we introduce a Spatial Transformer module, that can be included into a standard neural 
network architecture to provide spatial transformation capabilities. The action of the spatial trans¬ 
former is conditioned on individual data samples, with the appropriate behaviour learnt during train¬ 
ing for the task in question (without extra supervision). Unlike pooling layers, where the receptive 
fields are fixed and local, the spatial transformer module is a dynamic mechanism that can actively 
spatially transform an image (or a feature map) by producing an appropriate transformation for each 
input sample. The transformation is then performed on the entire feature map (non-locally) and 
can include scaling, cropping, rotations, as well as non-rigid deformations. This allows networks 
which include spatial transformers to not only select regions of an image that are most relevant (at¬ 
tention), but also to transform those regions to a canonical, expected pose to simplify recognition in 
the following layers. Notably, spatial transformers can be trained with standard back-propagation, 
allowing for end-to-end training of the models they are injected in. 
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Figure 1: The result of using a spatial transformer as the 
first layer of a fully-connected network trained for distorted 
MNIST digit classification, (a) The input to the spatial trans¬ 
former network is an image of an MNIST digit that is dis¬ 
torted with random translation, scale, rotation, and clutter, (b) 
The localisation network of the spatial transformer predicts a 
transformation to apply to the input image, (c) The output 
of the spatial transformer, after applying the transformation, 
(d) The classification prediction produced by the subsequent 
fully-connected network on the output of the spatial trans¬ 
former. The spatial transformer network (a CNN including a 
spatial transformer module) is trained end-to-end with only 
class labels - no knowledge of the groundtruth transforma¬ 
tions is given to the system. 


Spatial transformers can be incorporated into CNNs to benefit multifarious tasks, for example: 
(i) image classification: suppose a CNN is trained to perform multi-way classification of images 
according to whether they contain a particular digit - where the position and size of the digit may 
vary significantly with each sample (and are uncorrelated with the class); a spatial transformer that 
crops out and scale-normalizes the appropriate region can simplify the subsequent classification 
task, and lead to superior classification performance, see Fig. (ii) co-localisation: given a set of 
images containing different instances of the same (but unknown) class, a spatial transformer can be 
used to localise them in each image; (iii) spatial attention: a spatial transformer can be used for 
tasks requiring an attention mechanism, such as in but is more fiexible and can be trained 

purely with backpropagation without reinforcement learning. A key benefit of using attention is that 
transformed (and so attended), lower resolution inputs can be used in favour of higher resolution 
raw inputs, resulting in increased computational efficiency. 

The rest of the paper is organised as follows: Sect, [^discusses some work related to our own, we 
introduce the formulation and implementation of the spatial transformer in Sect. and finally give 
the results of experiments in Sect. Additional experiments and implementation details are given 
in Appendix]^ 

2 Related Work 

In this section we discuss the prior work related to the paper, covering the central ideas of modelling 
transformations with neural networks iniiisiiMi, learning and analysing transformation-invariant 
representations E 0 m Eoi Ea [331, as well as attention and detection mechanisms for feature 
selection EcmiKiiiiTiiia. 

Early work by Hinton ca looked at assigning canonical frames of reference to object parts, a theme 
which recurred in CD where 2D affine transformations were modeled to create a generative model 
composed of transformed parts. The targets of the generative training scheme are the transformed 
input images, with the transformations between input images and targets given as an additional 
input to the network. The result is a generative model which can learn to generate transformed 
images of objects by composing parts. The notion of a composition of transformed parts is taken 
further by Tieleman |[36l . where learnt parts are explicitly affine-transformed, with the transform 
predicted by the network. Such generative capsule models are able to learn discriminative features 
for classification from transformation supervision. 

The invariance and equivariance of CNN representations to input image transformations are studied 
in 1221 by estimating the linear relationships between representations of the original and transformed 
images. Cohen & Welling i) analyse this behaviour in relation to symmetry groups, which is also 
exploited in the architecture proposed by Gens & Domingos Co), resulting in feature maps that are 
more invariant to symmetry groups. Other attempts to design transformation invariant representa¬ 
tions are scattering networks IH, and CNNs that construct filter banks of transformed filters 1^13^ . 
Stollenga et al. nn use a policy based on a network’s activations to gate the responses of the net¬ 
work’s filters for a subsequent forward pass of the same image and so can allow attention to specific 
features. In this work, we aim to achieve invariant representations by manipulating the data rather 
than the feature extractors, something that was done for clustering in 191 . 

Neural networks with selective attention manipulate the data by taking crops, and so are able to learn 
translation invariance. Work such as GiEa are trained with reinforcement learning to avoid the 


2 














spatial Transformer 


Figure 2: The architecture of a spatial transformer module. The input feature map U is passed to a localisation 
network which regresses the transformation parameters 6. The regul ar sp atial grid G over V is transformed to 
the sampling grid TeiG), which is applied to U as described in Sect. 3.3 producing the warped output feature 
map V. The combination of the localisation network and sampling mechanism defines a spatial transformer. 


need for a differentiable attention mechanism, while ca use a differentiable attention mechansim 
by utilising Gaussian kernels in a generative model. The work by Girshick et al. CD uses a region 
proposal algorithm as a form of attention, and (71 show that it is possible to regress salient regions 
with a CNN. The framework we present in this paper can be seen as a generalisation of differentiable 
attention to any spatial transformation. 

3 Spatial Transformers 

In this section we describe the formulation of a spatial transformer. This is a differentiable module 
which applies a spatial transformation to a feature map during a single forward pass, where the 
transformation is conditioned on the particular input, producing a single output feature map. For 
multi-channel inputs, the same warping is applied to each channel. For simplicity, in this section we 
consider single transforms and single outputs per transformer, however we can generalise to multiple 
transformations, as shown in experiments. 

The spatial transformer mechanism is split into three parts, shown in Fig.[^ In order of computation, 
first a localisation network (Sect. HD takes the input feature map, and through a number of hidden 
layers outputs the parameters of the spatial transformation that should be applied to the feature map 
- this gives a transformation conditional on the input. Then, the predicted transformation parameters 
are used to create a sampling grid, which is a set of points where the input map should be sampled to 
produce the transformed output. This is done by the grid generator, described in Sect. |3.2[ Finally, 
the feature map and the sampling grid are take n as i nputs to the sampler, producing the output map 
sampled from the input at the grid points (Sect. |3.3| ). 

The combination of these three components forms a spatial transformer and will now be described 
in more detail in the following sections. 

3.1 Localisation Network 

The localisation network takes the input feature map U G width W, height H and 

C channels and outputs 0, the parameters of the transformation Te to be applied to the feature map: 
0 = f\oc{U). The size of 0 can vary depending on the transformation type that is parameterised, 
e.g. for an affine transformation 0 is 6-dimensional as in ( p^ . 

The localisation network function /ioc() can take any form, such as a fully-connected network or 
a convolutional network, but should include a final regression layer to produce the transformation 
parameters 0. 

3.2 Parameterised Sampling Grid 

To perform a warping of the input feature map, each output pixel is computed by applying a sampling 
kernel centered at a particular location in the input feature map (this is described fully in the next 
section). By pixel we refer to an element of a generic feature map, not necessarily an image. In 
general, the output pixels are defined to lie on a regular grid G = of pixels Gi = {x\^y\), 
forming an output feature map V G xVF xc^ where H' and W' are the height and width of the 
grid, and G is the number of channels, which is the same in the input and output. 
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Figure 3: Two examples of applying the parameterised sampling grid to an image U producing the output V. 
(a) The sampling grid is the regular grid G = Ti(G), where / is the identity transformation parameters, (b) 
The sampling grid is the result of warping the regular grid with an affine transformation TeiG). 

For clarity of exposition, assume for the moment that 7^ is a 2D affine transformation kg. We will 
discuss other transformations below. In this affine case, the pointwise transformation is 


= TeiGi) = ke | y\ 
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Vi 
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( 1 ) 


where yj) are the target coordinates of the regular grid in the output feature map, are 

the source coordinates in the input feature map that define the sample points, and Aq is the affine 
transformation matrix. We use height and width normalised coordinates, such that — 1 < < 1 

when within the spatial bounds of the output, and — 1 < xf^yf < 1 when within the spatial bounds 
of the input (and similarly for the y coordinates). The source/target transformation and sampling is 
equivalent to the standard texture mapping and coordinates used in graphics 10. 


The transform defined in ( p^ allows cropping, translation, rotation, scale, and skew to be applied 
to the input feature map, and requires only 6 parameters (the 6 elements of A^) to be produced by 
the localisation network. It allows cropping because if the transformation is a contraction (i.e. the 
determinant of the left 2x2 sub-matrix has magnitude less than unity) then the mapped regular grid 
will lie in a parallelogram of area less than the range of The effect of this transformation on 

the grid compared to the identity transform is shown in Fig.[^ 


The class of transformations Te may be more constrained, such as that used for attention 
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allowing cropping, translation, and isotropic scaling by varying s, tx, and ty. The transformation 
Te can also be more general, such as a plane projective transformation with 8 parameters, piece- 
wise affine, or a thin plate spline. Indeed, the transformation can have any parameterised form, 
provided that it is differentiable with respect to the parameters - this crucially allows gradients to be 
backpropagated through from the sample points Te{Gi) to the localisation network output 0. If the 
transformation is parameterised in a structured, low-dimensional way, this reduces the complexity 
of the task assigned to the localisation network. For instance, a generic class of structured and dif¬ 
ferentiable transformations, which is a superset of attention, affine, projective, and thin plate spline 
transformations, is 7^ = MeB, where 5 is a target grid representation (e.g. in B is the regu¬ 
lar grid G in homogeneous coordinates), and Me is a matrix parameterised by 0. In this case it is 
possible to not only learn how to predict 0 for a sample, but also to learn B for the task at hand. 


3.3 Differentiable Image Sampling 


To perform a spatial transformation of the input feature map, a sampler must take the set of sampling 
points Te{G), along with the input feature map U and produce the sampled output feature map V. 

Each (xf^yf) coordinate in Te{G) defines the spatial location in the input where a sampling kernel 
is applied to get the value at a particular pixel in the output V. This can be written as 

H W 

= E E - n; ^y) Vi G [1... H'W] Vc G [1... C] (3) 
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where and are the parameters of a generic sampling kernel k{) which defines the image 
interpolation (e.g. bilinear), is the value at location (n, m) in channel c of the input, and Vf 
is the output value for pixel i at location in channel c. Note that the sampling is done 

identically for each channel of the input, so every channel is transformed in an identical way (this 
preserves spatial consistency between channels). 

In theory, any sampling kernel can be used, as long as (sub-)gradients can be defined with respect to 
xl and yl. For example, using the integer sampling kernel reduces ^ to 

H W 

+ 0-5J - m)(5([y| + 0.5J - n) (4) 

n m 

where [x + 0.5j rounds x to the nearest integer and (^() is the Kronecker delta function. This 
sampling kernel equates to just copying the value at the nearest pixel to (x|, yf ) to the output location 
{xj, y\). Alternatively, a bilinear sampling kernel can be used, giving 

H W 

max(0,1 — |x| — m\) max(0,1 — \yl — n|) (5) 


To allow backpropagation of the loss through this sampling mechanism we can define the gradients 
with respect to U and G. For bilinear sampling <0 the partial derivatives are 


dVf 


H W 

EE max(0,1 — |x| — m\) max(0 ,1 — \y^ — n\) 

n m 


( 6 ) 


dxf 


H W 

E E max(0,1 - |y| - n|) 


n m 



if |m — x|| > 1 
if m > x| 
if m < x| 


(7) 


and similarly to (|7]) for . 

This gives us a (sub-)differentiable sampling mechanism, allowing loss gradients to flow back not 
only to the input feature map but also to the sampling grid coordinates d^, and therefore back 

dx^ doc^ 

to the transformation parameters 0 and localisation network since and can be easily derived 
from for example. Due to discontinuities in the sampling factions, sub-gradients must be used. 
This sampling mechanism can be implemented very efficiently on GPU, by ignoring the sum over 
all input locations and instead just looking at the kernel support region for each output pixel. 

3.4 Spatial Transformer Networks 


The combination of the localisation network, grid generator, and sampler form a spatial transformer 
(Fig-ig. This is a self-contained module which can be dropped into a CNN architecture at any point, 
and in any number, giving rise to spatial transformer networks. This module is computationally very 
fast and does not impair the training speed, causing very little time overhead when used naively, and 
even speedups in attentive models due to subsequent downsampling that can be applied to the output 
of the transformer. 


Placing spatial transformers within a CNN allows the network to learn how to actively transform 
the feature maps to help minimise the overall cost function of the network during training. The 
knowledge of how to transform each training sample is compressed and cached in the weights of 
the localisation network (and also the weights of the layers previous to a spatial transformer) during 
training. For some tasks, it may also be useful to feed the output of the localisation network, 0, 
forward to the rest of the network, as it explicitly encodes the transformation, and hence the pose, of 
a region or object. 

It is also possible to use spatial transformers to downsample or oversample a feature map, as one can 
define the output dimensions H' and W' to be different to the input dimensions H and W. However, 
with sampling kernels with a fixed, small spatial support (such as the bilinear kernel), downsampling 
with a spatial transformer can cause aliasing effects. 
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MNIST Distortion 

Model 

R 

RTS 

P 

E 

FCN 

2.1 

5.2 

3.1 

3.2 

CNN 

1.2 

0.8 

1.5 

1.4 

Aff 

1.2 

0.8 

1.5 

2.7 

ST-FCN Proj 

1.3 

0.9 

1.4 

2.6 

TPS 

1.1 

0.8 

1.4 

2.4 

Aff 

0.7 

0.5 

0.8 

1.2 

ST-CNN Proj 

0.8 

0.6 

0.8 

1.3 

TPS 

0.7 

0.5 

0.8 

1.1 




Table 1: Left: The percentage errors for different models on different distorted MNIST datasets. The different 
distorted MNIST datasets we test are TC: translated and cluttered, R: rotated, RTS: rotated, translated, and 
scaled, P: projective distortion, E: elastic distortion. All the models used for each experiment have the same 
number of parameters, and same base structure for all experiments. Right: Some example test images where 
a spatial transformer network correctly classifies the digit but a CNN fails, (a) The inputs to the networks, (b) 
The transformations predicted by the spatial transformers, visualised by the grid Tq (G). (c) The outputs of the 
spatial transformers. E and RTS examples use thin plate spline spatial transformers (ST-CNN TPS), while R 
examples use affine spatial transformers (ST-CNN Aff) with the angles of the affine transformations given. Eor 
videos showing animations of these experiments and more see https : / /goo . gl/qdEhUu 

Finally, it is possible to have multiple spatial transformers in a CNN. Placing multiple spatial trans¬ 
formers at increasing depths of a network allow transformations of increasingly abstract representa¬ 
tions, and also gives the localisation networks potentially more informative representations to base 
the predicted transformation parameters on. One can also use multiple spatial transformers in paral¬ 
lel - this can be useful if there are multiple objects or parts of interest in a feature map that should be 
focussed on individually. A limitation of this architecture in a purely feed-forward network is that 
the number of parallel spatial transformers limits the number of objects that the network can model. 


4 Experiments 


In this section we expl ore the use of spatial transformer networks on a number of supervised learn¬ 
ing tasks. In Sect. |4.1| we begin with experiments on distorted versions of the MNIST handwriting 
dataset, showing the ability of spatial transformers to improve classification performance through 
actively transforming the input images. In Sect. 4.2 we test spatial transformer networks on a chal¬ 
lenging real-world dataset, Street View House Numbers (251, foi* number recognition, showing state- 
of-the-art results using multiple spatial transformers embedded in the convolutional stack of a CNN. 
Finally, in Sect. |4.3[ we investigate the use of multiple parallel spatial transformers for fine-grained 
classification, showing state-of-the-art performance on CUB-200-2011 birds dataset (3^ by dis¬ 
covering object parts and learning to attend to them. Further experiments of MNIST addition and 
co-localisation can be found in Appendix [A[ 

4.1 Distorted MNIST 


In this section we use the MNIST handwriting dataset as a testbed for exploring the range of trans¬ 
formations to which a network can learn invariance to by using a spatial transformer. 

We begin with experiments where we train different neural network models to classify MNIST data 
that has been distorted in various ways: rotation (R), rotation, scale and translation (RTS), projective 
transformation (P), and elastic warping (E) - note that elastic warping is destructive and can not be 
inverted in some cases. The full details of the distortions used to generate this data are given in 
Appendix]^ We train baseline fully-connected (FCN) and convolutional (CNN) neural networks, 
as well as networks with spatial transformers acting on the input before the classification network 
(ST-FCN and ST-CNN). The spatial transformer networks all use bilinear sampling, but variants use 
different transformation functions: an affine transformation (Aff), projective transformation (Proj), 
and a 16-point thin plate spline transformation (TPS) O. The CNN models include two max-pooling 
layers. All networks have approximately the same number of parameters, are trained with identical 
optimisation schemes (backpropagation, SGD, scheduled learning rate decrease, with a multinomial 
cross entropy loss), and all with three weight layers in the classification network. 

The results of these experiments are shown in Table(left). Looking at any particular type of dis¬ 
tortion of the data, it is clear that a spatial transformer enabled network outperforms its counterpart 
base network. For the case of rotation, translation, and scale distortion (RTS), the ST-CNN achieves 
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Size 

Model 

64px 

128px 

Maxout CNNll3i 

4.0 

- 

CNN (ours) 

4.0 

5.6 

DRAM* C] 

3.9 

4.5 

ST-CNN 

Multi 

3.7 

3.6 

3.9 

3.9 



ST-^ conv -►ST-*- conv -^ST^ conv -► 

. 


ST-^-- 



2 

6 
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(b) 
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Table 2: Left: The sequence error for SVHN multi-digit recognition on crops of 64 x 64 pixels (64px), and 
inflated crops of 128 x 128 (128px) which include more background. *The best reported result from [T] uses 
model averaging and Monte Carlo averaging, whereas the results from other models are from a single forward 
pass of a single model. Right: (a) The schematic of the ST-CNN Multi model. The transformations applied by 
each spatial transformer (ST) is applied to the convolutional feature map produced by the previous layer, (b) 
The result of multiplying out the affine transformations predicted by the four spatial transformers in ST-CNN 
Multi, visualised on the input image. 


0.5% and 0.6% depending on the class of transform used for Te, whereas a CNN, with two max¬ 
pooling layers to provide spatial invariance, achieves 0.8% error. This is in fact the same error that 
the ST-FCN achieves, which is without a single convolution or max-pooling layer in its network, 
showing that using a spatial transformer is an alternative way to achieve spatial invariance. ST-CNN 
models consistently perform better than ST-FCN models due to max-pooling layers in ST-CNN pro¬ 
viding even more spatial invariance, and convolutional layers better modelling local structure. We 
also test our models in a noisy environment, on 60 x 60 images with translated MNIST digits and 
background clutter (see Fig. [T] third row for an example): an FCN gets 13.2% error, a CNN gets 
3.5% error, while an ST-FCN gets 2.0% error and an ST-CNN gets 1.7% error. 

Looking at the results between different classes of transformation, the thin plate spline transfor¬ 
mation (TPS) is the most powerful, being able to reduce error on elastically deformed digits by 
reshaping the input into a prototype instance of the digit, reducing the complexity of the task for the 
classification network, and does not over fit on simpler data e.g. R. Interestingly, the transformation 
of inputs for all ST models leads to a “standard” upright posed digit - this is the mean pose found 
in the training data. In Table (right), we show the transformations performed for some test cases 
where a CNN is unable to correctly classify the digit, but a spatial transformer network can. Further 
test examples are visualised in an animation here https : / / goo . gl/qdEhUu 

4.2 Street View House Numbers 

We now test our spatial transformer networks on a challenging real-world dataset, Street View House 
Numbers (SVHN) 1^ . This dataset contains around 200k real world images of house numbers, with 
the task to recognise the sequence of numbers in each image. There are between 1 and 5 digits in 
each image, with a large variability in scale and spatial arrangement. 

We follow the experimental setup as in Elia, where the data is preprocessed by taking 64 x 64 
crops around each digit sequence. We also use an additional more loosely 128 x 128 cropped dataset 
as in ID. We train a baseline character sequence CNN model with 11 hidden layers leading to five 
independent softmax classifiers, each one predicting the digit at a particular position in the sequence. 
This is the character sequence model used in (191, where each classifier includes a null-character 
output to model variable length sequences. This model matches the results obtained in m. 

We extend this baseline CNN to include a spatial transformer immediately following the input (ST- 
CNN Single), where the localisation network is a four-layer CNN. We also define another extension 
where before each of the first four convolutional layers of the baseline CNN, we insert a spatial 
transformer (ST-CNN Multi), where the localisation networks are all two layer fully connected net¬ 
works with 32 units per layer. In the ST-CNN Multi model, the spatial transformer before the first 
convolutional layer acts on the input image as with the previous experiments, however the subse¬ 
quent spatial transformers deeper in the network act on the convolutional feature maps, predicting a 
transformation from them and transforming these feature maps (this is visualised in Table (right) 
(a)). This allows deeper spatial transformers to predict a transformation based on richer features 
rather than the raw image. All networks are trained from scratch with SGD and dropout ITtI. with 
randomly initialised weights, except for the regression layers of spatial transformers which are ini¬ 
tialised to predict the identity transform. Affine transformations and bilinear sampling kernels are 
used for all spatial transformer networks in these experiments. 
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Model 


Cimpoi ’15 151 

66.7 

Zhang ’14 l40l 

74.9 

Branson ’14 O 

75.7 

Lin’15 

80.9 

Simon ’15 

81.0 

CNN (ours) 224px 

82.3 

2 X ST-CNN 224px 

83.1 

2 X ST-CNN 448px 

83.9 

4 X ST-CNN 448px 

84.1 



Table 3: Left: The accuracy on CUB-200-2011 bird classification dataset. Spatial transformer networks with 
two spatial transformers (2 x ST-CNN) and four spatial transformers (4 x ST-CNN) in parallel achieve higher 
accuracy. 448px resolution images can be used with the ST-CNN without an increase in computational cost 
due to downsampling to 224px after the transformers. Right: The transformation predicted by the spatial 
transformers of 2 x ST-CNN (top row) and 4 x ST-CNN (bottom row) on the input image. Notably for the 
2 X ST-CNN, one of the transformers (shown in red) learns to detect heads, while the other (shown in green) 
detects the body, and similarly for the 4 x ST-CNN. 

The results of this experiment are shown in Table (left) - the spatial transformer models obtain 
state-of-the-art results, reaching 3.6% error on 64 x 64 images compared to previous state-of-the-art 
of 3.9% error. Interestingly on 128 x 128 images, while other methods degrade in performance, 
an ST-CNN achieves 3.9% error while the previous state of the art at 4.5% error is with a recurrent 
attention model that uses an ensemble of models with Monte Carlo averaging - in contrast the ST- 
CNN models require only a single forward pass of a single model. This accuracy is achieved due to 
the fact that the spatial transformers crop and rescale the parts of the feature maps that correspond 
to the digit, focussing resolution and network capacity only on these areas (see Table (right) (b) 
for some examples). In terms of computation speed, the ST-CNN Multi model is only 6% slower 
(forward and backward pass) than the CNN. 


4.3 Fine-Grained Classification 


In this section, we use a spatial transformer network with multiple transformers in parallel to perform 
fine-grained bird classification. We evaluate our models on the CUB-200-2011 birds dataset |[38l, 
containing 6k training images and 5.8k test images, covering 200 species of birds. The birds appear 
at a range of scales and orientations, are not tightly cropped, and require detailed texture and shape 
analysis to distinguish. In our experiments, we only use image class labels for training. 

We consider a strong baseline CNN model - an Inception architecture with batch normalisation ca 
pre-trained on ImageNet and fine-tuned on CUB - which by itself achieves the state-of-the- 
art accuracy of 82.3% (previous best result is 81.0% Eol). We then train a spatial transformer 
network, ST-CNN, which contains 2 or 4 parallel spatial transformers, parameterised for attention 
and acting on the input image. Discriminative image parts, captured by the transformers, are passed 
to the part description sub-nets (each of which is also initialised by Inception). The resulting part 
representations are concatenated and classified with a single softmax layer. The whole architecture 
is trained on image class labels end-to-end with backpropagation (full details in Appendix [A|. 

The results are shown in Table(left). The ST-CNN achieves an accuracy of 84.1%, outperforming 
the baseline by 1.8%. It should be noted that there is a small (22/5794) overlap between the Ima¬ 
geNet training set and CUB-200-2011 test se|^- removing these images from the test set results in 
84.0% accuracy with the same ST-CNN. In the visualisations of the transforms predicted by 2 x ST- 
CNN (Table (right)) one can see interesting behaviour has been learnt: one spatial transformer 
(red) has learnt to become a head detector, while the other (green) fixates on the central part of the 
body of a bird. The resulting output from the spatial transformers for the classification network is 
a somewhat pose-normalised representation of a bird. While previous work such as ID explicitly 
define parts of the bird, training separate detectors for these parts with supplied keypoint training 
data, the ST-CNN is able to discover and learn part detectors in a data-driven manner without any 
additional supervision. In addition, the use of spatial transformers allows us to use 448px resolution 
input images without any impact in performance, as the output of the transformed 448px images are 
downsampled to 224px before being processed. 


^Thanks to the eagle-eyed Hugo Larochelle and Yin Zheng for spotting the birds nested in both the ImageNet 
training set and CUB test set. 
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5 Conclusion 


In this paper we introduced a new self-contained module for neural networks - the spatial trans¬ 
former. This module can be dropped into a network and perform explicit spatial transformations 
of features, opening up new ways for neural networks to model data, and is learnt in an end-to- 
end fashion, without making any changes to the loss function. While CNNs provide an incredibly 
strong baseline, we see gains in accuracy using spatial transformers across multiple tasks, result¬ 
ing in state-of-the-art performance. Furthermore, the regressed transformation parameters from the 
spatial transformer are available as an output and could be used for subsequent tasks. While we 
only explore feed-forward networks in this work, early experiments show spatial transformers to be 
powerful in recurrent models, and useful for tasks requiring the disentangling of object reference 
frames, as well as easily extendable to 3D transformations (see Appendix ]A. 3| ). 
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A Appendix 

In this section we present the results of two further experiments - that of MNIST addition showing 
spatial transformers acting on multiple objects in Sect. |A.l[ and co-localisation in Sect. |A^ showing 
the applicat ion t o semi-supervised scenarios. In addition, we give an example o f the extensi on to 
3D in Sect. |A.3[ We also expand upon the details of the experiments from Sect. |4.1| in Sect. |A.4| 
Sect. |4.2| in Sect. |A.5| and Sect. |4.3| in Sect. |A.6| 

A.l MNIST Addition 

In this section we demonstrate another use case for multiple spatial transformers in parallel: to 
model multiple objects. We define an MNIST addition task, where the network must output the sum 
of the two digits given in the input. Each digit is presented in a separate 42 x 42 input channel 
(giving 2-channel inputs), but each digit is transformed independently, with random rotation, scale, 
and translation (RTS). 

We train fully connected (FCN), convolutional (CNN) and single spatial transformer fully connected 
(ST-FCN) networks, as well as spatial transformer fully connected networks with two parallel spa¬ 
tial transformers (2 x ST-FCN) acting on the input image, each one taking both channels as input and 
transforming both channels. The two 2-channel outputs of the two spatial transformers are concate¬ 
nated into a 4-channel feature map for the subsequent FCN. As in Sect. |4.1[ all networks have the 
same number of parameters, and are all trained with SGD to minimise the multinomial cross entropy 
loss for 19 classes (the possible addition results 0-18). 

The results are given in Table (left). Due to the complexity of this task, the FCN reaches a 
minimum error of 47.7%, however a CNN with max-pooling layers is far more accurate with 14.7% 
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Concatenated 


Model 

RTS 

FCN 

47.7 

CNN 

14.7 

Aff 

22.6 

ST-FCN Proj 

18.5 

TPS 

19.1 

Aff 

9.0 

2 X ST-FCN Proj 

5.9 

TPS 

5.8 



Table 4: Left: The percentage error for the two digit MNIST addition task, where each digit is transformed 
independently in separate channels, trained by supplying only the label of the sum of the two digits. The 
use of two spatial transformers in parallel, 2 x ST-FCN, allows the fully-connected neural network to become 
invariant to the transformations of each digit, giving the lowest error. All the models used for each column 
have approximately the same number of parameters. Right: A test example showing the learnt behaviour of 
each spatial transformer (using a thin plate spline (TPS) transformation). The 2-channel input (the blue bar 
denotes separation between channels) is fed to two independent spatial transformers, STl and ST2, each of 
which operate on both channels. The outputs of STl and ST2 and concatenated and used as a 4-channel input 
to a fully connected network (FCN) which predicts the addition of the two original digits. During training, the 
two spatial transformers co-adapt to focus on a single channel each. 



MNIST Distortion 

Class 

T 

TC 

0 

100 

81 

1 

100 

82 

2 

100 

88 

3 

100 

75 

4 

100 

94 

5 

100 

84 

6 

100 

93 

7 

100 

85 

8 

100 

89 

9 

100 

87 


In 


In 



<ll) 

) small 
distance 

e{ll) 

5 large 

distanee 


Table 5: Left: The percent of correctly co-localised digits for different MNIST digit classes, for just translated 
digits (T), and for translated digits with clutter added (TC). Right: The optimisation architecture. We use a 
hinge loss to enforce the distance between the two outputs of the spatial transformer (ST) to be less than the 
distance to a random crop, hoping to encourage the spatial transformer to localise the common objects. 


error. Adding a single spatial transformer improves the capability of an FCN by focussing on a 
single region of the input containing both digits, reaching 18.5% error. However, by using two 
spatial transformers, each transformer can learn to focus on transforming the digit in a single channel 
(though receiving both channels as input), visualised in Table(right). The transformers co-adapt, 
producing stable representations of the two digits in two of the four output channels of the spatial 
transformers. This allows the 2 xST-FCN model to achieve 5.8% error, far exceeding that of other 
models. 

A.2 Co-localisation 

In this experiment, we explore the use of spatial transformers in a semi-supervised scenario - co¬ 
localisation. The co-localisation task is as follows: given a set of images that are assumed to contain 
instances of a common but unknown object class, localise (with a bounding box) the common object. 
Neither the object class labels, nor the object location ground truth is used for optimisation, only the 
set of images. 

To achieve this, we adopt the supervision that the distance between the image crop corresponding to 
two correctly localised objects is smaller than to a randomly sampled image crop, in some embed¬ 
ding space. For a dataset X = {In} of N images, this translates to a triplet loss, where we minimise 
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Optimisation 



Figure 4: A look at the optimisation dynamics for co-localisation. Here we show the localisation predicted 
by the spatial transformer for three of the 100 dataset images after the SGD step labelled below. By SGD 
step 180 the model has process has correctly localised the three digits. A full animation is shown in the video 

https://goo.gl/qdEhUu 

the hinge loss 


N M 

Y, Y max(0, ||e(Jt) - e{ll)\\l - \\e{lZ) - e(/r)||^ +a) (8) 

n m^n 


where is the image crop of In corresponding to the localised object, is a randomly sampled 
patch from e() is an encoding function and a is a margin. We can use a spatial transformer to 
act as the localises such that = Te{In) where 0 = f\oc{In)^ interpreting the parameters of the 
transformation 0 as the bounding box of the object. We can minimise this with stochastic gradient 
descent, randomly sampling image pairs (n, m). 


We perform co-localisation on translated (T), and also translated and cluttered (TC) MNIST images. 
Each image, a 28 x 28 pixel MNIST digit, is placed in a uniform random location in a 84 x 84 
black background image. For the cluttered dataset, we also then add 16 random 6x6 crops sam¬ 
pled from the original MNIST training dataset, creating distractors. For a particular co-localisation 
optimisation, we pick a digit class and generate 100 distorted image samples as the dataset for the 
experiment. We use a margin a = 1, and for the encoding function e() we use the CNN trained for 
digit classification from Sect. 4.1 concatenating the three layers of activations (two hidden layers 
and the classification layer without softmax) to form a feature descriptor. We use a spatial trans¬ 
former parameterised for attention (scale and translation) where the localisation network is a 100k 
parameter CNN consisting of a convolutional layer with eight 9x9 filters and a 4 pixel stride, fol¬ 
lowed by 2 X 2 max pooling with stride 2 and then two 8-unit fully-connected layers before the final 
3-unit fully-connected layer. 


The results are shown in TableWe measure a digit to be correctly localised if the overlap (area of 
intersection divided by area of union) between the predicted bounding box and groundtruth bounding 
box is greater than 0.5. Our co-localisation framework is able to perfectly localise MNIST digits 
without any clutter with 100% accuracy, and correctly localises between 75-93% of digits when 
there is clutter in the images. An example of the optimisation process on a subset of the dataset for 
“8” is shown in Fig. This is surprisingly good performance for what is a simple loss function 
derived from simple intuition, and hints at potential further applications in tracking problems. 


A.3 Higher Dimensional Transformers 


The framework described in this paper is not limited to 2D transformations and can be easily ex¬ 
tended to higher dimensions. To demonstrate this, we give the example of a spatial transformer 
capable of performing 3D affine transformations. 
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3D transformation applied 



Figure 5: The behaviour of a trained 3D MNIST classifier on a test example. The 3D voxel input contains 
a random MNIST digit which has been extruded and randomly placed inside a 60 x 60 x 60 volume. A 3D 
spatial transformer performs a transformation of the input, producing an output volume whose depth is then 
fiattened. This creates a 2D projection of the 3D space, which the subsequent layers of the network are able to 
classify. The whole network is trained end-to-end with just classification labels. 


We extended the differentiable image sampling of Sect. 3.3 to perform 3D bilinear sampling. The 
3D equivalent of ^ becomes 


H W D 

max(0,1 - \xt - m\) inax(0,1 - |y| - n|) max(0,1 - \zt - ^|) (9) 

n m I 


for the 3D input U G and output V G r^'xw'xd'xC^ W', and B' are the 

height, widt h an d depth of the grid, and C is the number of channels. Similarly to the 2D sampling 
grid in Sect. 3.2 the source coordinates that define the sampling points, (x|, zf) can be generated 


by the transformation of a regular 3D grid G = of voxels Gi = For a 3D affine 

transformation this is 
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The 3D spatial transformer can be used just like its 2D counterpart, being dropped into neural net¬ 
works to provide a way to warp data in 3D space, where the third dimension could be space or 
time. 


Another interesting way to use the 3D transformer is to fiatten the 3D output across one dimension, 
creating a 2D projection of the 3D space, e.g. ynmi ^ ^ ^h'xw'xC^ 

allows the original 3D data to be intelligently projected to 2D, greatly reducing the dimensionality 
and complexity of the subsequent processing. We demonstrated this on the task of 3D object classi¬ 
fication on a dataset of 3D, extruded MNIST digits. The task is to take a 3D voxel input of a digit 
which has been randomly translated and rotated in 3D space, and output the class of the digit. The 
resulting 3D spatial transformer network learns to create a 2D projection of the 3D space where the 
digit is centered in the resulting 2D image, making it easy for the remaining layers to classify. An 
example is shown in Fig.[^ 

A.4 Distorted MNIST Details 


In this section we expand upon the details of the distorted MNIST experiments in Sect. 4.1 


Data. The rotated dataset (R) was generated from rotating MNIST training digits with a random 
rotation sampled uniformly between —90° and +90°. The rotated, translated, and scaled dataset 
(RTS) was generated by randomly rotating an MNIST digit by +45° and —45°, randomly scaling the 
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digit by a factor of between 0.7 and 1.2, and placing the digit in a random location in a 42 x 42 image, 
all with uniform distributions. The projected dataset (P) was generated by scaling a digit randomly 
between 0.75 and 1.0, and stretching each comer of an MNIST digit by an amount sampled from a 
normal distribution with zero mean and 5 pixel standard deviation. The elasticly distorted dataset 
(E) was generated by scaling a digit randomly between 0.75 and 1.0, and then randomly peturbing 
16 control points of a thin plate spline arranged in a regular grid on the image by an amount sampled 
from a normal distribution with zero mean and 1.5 pixel standard deviation. The translated and 
cluttered dataset (TC) is generated by placing an MNIST digit in a random location in a 60 x 60 
black canvas, and then inserting six randomly sampled 6x6 patches of other digit images into 
random locations in the image. 

Networks. All networks use rectified linear non-linearities and softmax classifiers. All FCN net¬ 
works have two hidden fully connected layers followed by a classification layer. All CNN networks 
have a 9 X 9 convolutional layer (stride 1, no padding), a 2 x 2 max-pooling layer with stride 2, a 
subsequent 7x7 convolutional layer (stride 1, no padding), and another 2x2 max-pooling layer 
with stride 2 before the final classfication layer. All spatial transformer (ST) enabled networks place 
the ST modules at the beginning of the network, and have three hidden layers in their localisation 
networks with 32 unit fully connected layers for ST-FCN networks and two 20-filter 5x5 convo¬ 
lutional layers (stride 1, no padding) acting on a 2x downsampled input, with 2x2 max-pooling 
between convolutional layers, and a 20 unit fully connected layer following the convolutional layers. 
Spatial transformer networks for TC and RTS datasets have average pooling after the spatial trans¬ 
former to downsample the output of the transformer by a factor of 2 for the classification network. 
The exact number of units in FCN and CNN based classification models varies so as to always en¬ 
sure that all networks for a particular experiment contain the same number of learnable parameters 
(around 400k). This means that spatial transformer networks generally have less parameters in the 
classification networks due to the need for parameters in the localisation networks. The FCNs have 
between 128 and 256 units per layer, and the CNNs have between 32 and 64 filters per layer. 

Training. All networks were trained with SGD for 150k iterations, the same hyperparameters (256 
batch size, 0.01 base learning rate, no weight decay, no dropout), and same learning rate schedule 
(learning rate reduced by a factor of ten every 50k iterations). We initialise the network weights ran¬ 
domly, except for the final regression layer of localisation networks which are initialised to regress 
the identity transform (zero weights, identity transform bias). We perform three complete training 
runs for all models with different random seeds and report average accuracy. 

A.5 Street View House Numbers Details 


For the SVHN experiments in Sect. |4.2[ we follow (TJ [T3J and select hyperparameters from a vali¬ 
dation set of 5k images from the training set. All networks are trained for 400k iterations with SGD 
(128 batch size), using a base learning rate of 0.01 decreased by a factor of ten every 80k iterations, 
weight decay set to 0.0005, and dropout at 0.5 for all layers except the first convolutional layer and 
localisation networks. The learning rate for localisation networks of spatial transformer networks 
was set to a tenth of the base learning rate. 


We adopt the notation that corwlN, w,s,p] denotes a convolutional layer with N filters of size 
w X w, with stride s and p pixel padding, fc[A^] is a fully connected layer with N units, and 
max[s] is a 5 X 5 max-pooling layer with stride s. The CNN model is: conv[48,5,l,2]-max[2]- 
conv[64,5,l,2]-conv[128,5,l,2]-max[2]-conv[160,5,l,2]-conv[192,5,l,2]-max[2]-conv[192,5,l,2]- 
conv[192,5,l,2]-max[2]-conv[192,5,l,2]-fc[3072]-fc[3072]-fc[3072], with rectified linear units 
following each weight layer, followed by five parallel fc[ll] and softmax layers for classification 
(similar to that in |[T9l ). The ST-CNN Single has a single spatial transformer (ST) before the 
first convolutional layer of the CNN model - the ST’s localisation network architecture is as 
follows: conv[32,5,l,2]-max[2]-conv[32,5,l,2]-fc[32]-fc[32]. The ST-CNN Multi has four spatial 
transformers, one before each of the first four convolutional layers of the CNN model, and each 
with a simple fc[32]-fc[32] localisation network. 

We initialise the network weights randomly, except for the final regression layer of localisation 
networks which are initialised to regress the identity transform (zero weights, identity transform 
bias). We performed two full training runs with different random seeds and report the average 
accuracy obtained by a single model. 
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224px 



224px 


Figure 6: The architecture of the 2 x ST-CNN 448px used for bird classification. A single localisation network 
/loc predicts two transformation parameters 0i and O 2 , with the subsequent transforms Te-i and 7^2 applied to 
the original input image. 

A.6 Fine Grained Classification Details 

In this section we describe our fine-grained image classification architecture in more detail. For 
this task, we utilise the spatial transformers as a differentiable attention mechanism, where each 
transformer is expected to automatically learn to focus on discriminative object parts. Namely, each 
transformer predicts the location (x,y) of the attention window, while the scale is fixed to 50% of the 
image size. The transformers sample 224 x 224 crops from the input image, each of which is then 
described each by its own CNN stream, thus forming a multi-stream architecture (shown in Fig.[^. 
The outputs of the streams are 1024-D crop descriptors, which are concatenated and classified with 
a 200-way softmax classifier. 

As the main building block of our network, we utilise the state-of-the-art Inception architecture with 
batch normalisation d, pre-trained on the ImageNet Challenge (ILSVRC) dataset. Our model 
achieves 27.1% top-1 error on the ILSVRC validation set using a single image crop (we only trained 
on single-scale images, resized so that the smallest side is 256). The crop description networks 
employ the Inception architecture with the last layer (1000-way ILSVRC classifier) removed, so 
that the output is a 1024-D descriptor. 

The localisation network is shared across all the transformers, and was derived from Inception in 
the following way. Apart from the ILSVRC classification layer, we also removed the last pooling 
layer to preserve the spatial information. The output of this truncated Inception net has 7x7 spatial 
resolution and 1024 feature channels. On top of it, we added three weight layers to predict the 
transformations: (i) 1 x 1 convolutional layer to reduce the number of feature channels from 1024 
to 128; (ii) fully-connected layer with 128-D output; (iii) fully-connected layer with 2A^-D output, 
where N is the number of transformers (we experimented with N = 2 and N = 4). 

We note that we did not strive to optimise the architecture in terms of the number of parameters 
and the computation time. Our aim was to investigate whether spatial transformer networks are 
able to automatically discover meaningful object parts when trained just on image labels, which we 
confirmed both quantitatively and qualitatively (Sect. |4.3| . 

The model was trained for 30k iterations with SGD (batch size 256) with an initial learning rate 
of 0.1, reduced by a factor of 10 after 10k, 20k, and 25k iterations. For stability, the localisation 
network’s learning rate is the base learning rate multiplied by 10“^. Weight decay was set at 10“ 5 
and dropout of 0.7 was used before the 200-way classification layer. 

We evaluated two input images sizes for the spatial transformers: 224 x 224 and 448 x 448. In the 
latter case, we added a fixed 2 x downscaling layer before the localisation net, so that its input is 
still 224 X 224. The difference between the two settings lies in the size of the image from which 
sampling is performed (224 vs 448), with 448 better suited for sampling small-scale crops. The 
output of the transformers are 224 x 224 crops in both cases (so that they are compatible with crop 
description Inception nets). When training, we utilised conventional augmentation in the form of 
random sampling (224 x 224 from 256 x S and 448 x 448 from 512 x S' where S is the largest 
image side) and horizontal flipping. The localisation net was initialised to tile the image plane with 
the spatial transformer crops. 

We also experimented with more complex transformations (location and scale, as well as affine), but 
observed similar results. This can be attributed to the very small size of the training set (6k images, 
200 classes), and we noticed severe over-fitting in all training scenarios. The hyper-parameters were 
estimated by cross-validation on the training set. 
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